Generation of entangled light via dynamical Casimir effect 
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This paper addresses the excitation of vacuum fluctuations of the electromagnetic field through 
periodic modulations of a refractive index and the possibility of using entanglement as a distinctive 
marker of the quantum nature of the phenomenon. It introduces a lossy environment and analyses its 
implications on the possibility of generating such an effect and measuring entanglement, concluding 
that it is not entirely destructive when the produced particles share the same environment. 
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Introduction. The dynamical Casimir effect (DCE) 
consists basically in applying time varying boundary con- 
ditions or refractive index modulations to the electromag- 
netic vacuum in order to excite real particles out of the 
zero-point fluctuations (ZPF) [IHl]. Even though this is 
probably the most direct method to probe and enhance 
ZPF, detecting it in optical cavities is not easy, albeit 
there being some promising attempts in that direction 
[SI |B] and observations in analogue experiments [7] . The 
obstacles have mainly been attributed to the difficulty in 
producing intense and fast enough optical perturbations 
[SI IHj and to technical aspects of the experiment (cavity 
losses, detector sensitivity, etc.) capable of measuring 
the scarce radiation produced [8l [9] . Finite temperature 
has two added consequences: it provides thermal photons 
which seed the DCE, thus masking the ZPF, and pro- 
duces decoherence and losses via the thermalisation of 
the optical cavity by the environment. In fact, previous 
works have shown that a seed thermal state greatly im- 
proves the photon production [TO]. However, as it will 
be demonstrated here, this amplification of thermal radi- 
ation is also predicted with a classical model, rendering 
the DCE at finite temperatures, in practice, a classical 
effect. 

The key aspect of zero temperature DCE is that, un- 
like thermal fluctuations, ZPF are coherent and the para- 
metric amplification produced by a time varying optical 
cavity preserves this coherence. For example, in an al- 
ternative DCE scheme, which considers a medium with a 
time varying refractive index in the absence of any spatial 
boundary, known as time refraction (TR) , the parametric 
amplification of ZPF generates pairs of counter propagat- 
ing photons which are maximally entangled and exit the 
medium separately [11]. For DCE in linear cavities the 
counter propagating photons reflect into each other, ren- 
dering them indistinguishable and thus impeding the use 
of entanglement as a marker of quantum character. How- 
ever, this difficulty can be circumvented using different 
geometries. 

Previous work in the problem considered the enhance- 
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ment of the output radiation due to thermal seeding [10] 
and the impact of losses [SUHldl] (phase damping, dissip- 
ation, Markovian baths, etc.) in linear optical cavities. 
At finite temperature (T > 0) the losses in the cavity 
do not force the cavity to the vacuum state but to the 
thermal state of the cavity which is in equilibrium with 
the environment. This provides constant thermal reseed- 
ing of the cavity (not just at the initial moment) masking 
even further any mark of ZPF. 

This letter introduces the idea of producing DCE in 
an optical ring resonator (like a Sagnac interferometer) 
filled with a dielectric subjected to a periodic modulation 
of the refractive index. This configuration possesses the 
following advantages: (i) the pairs of photons generated 
by DCE do not refiect into each other but propagate in 
opposite directions and can be extracted independently 
(see Fig. II]), remaining therefore distinguishable; (ii) the 
two photon beams propagate along the same optical path 
and experience the same decoherence and losses (they 
share the same bath) which allows for the existence of a 
decoherence free subspace (DFS) [13]; and (iii) the degree 
of entanglement between the two beams is a mark of ZPF 
even at finite temperature. 

The models for DCE in linear cavities predict that the 
number of photons only grows exponentially if the rate 
of their creation is faster than they decay (weak losses 
regime) l9l, but even then their growth rate is reduced 
by losses. For ring cavities, the DCE can support an 
exponential photon generation and entanglement, even 
in the strong losses regime. 

Classical vs. quantum model. The classical model 
for an optical cavity filled with a time dependent dielec- 
tric is basically the same for both linear and ring res- 
onators and can be derived directly from the macro- 
scopic Maxwell equations by imposing a time modula- 
tion of the refractive index n = n{t) — uq + Sn h{t), 
with h{t) E [—1, 1]. The resulting dynamical eq. for the 
amplitude of electric field E of an harmonic mode is 



E + 4rmE + (w^ + 2hn + 2n^)E = 



(1) 



and can be converted into the Hill equation [14] 

y = ~fy, (2) 




^i 



Figure 1: Linear and ring resonators 

Different cavity configurations. In tlie linear cavity the 

eigenmodes are standing waves and there is only one mode for 

each frequency. In the ring cavity two counterpropagating 

traveUing wave modes of the same frequency can coexist and be 

extracted independently. 



in eq. (IT]) by an operator and computing the effective 
Hamiltonian for each mode. The difference between Un- 
ear and ring resonators results from the mode structure 
of the E operator determined by the geometry of the cav- 
ity. The effects of losses and decoherence are modelled 
via the coupling to a Markovian bath according to the 
master eq. (similar to refs p,ll3j) 



P = -I 



{n 4- 1/2) [aU]P + pal 



+ n{al paj - p 



}(6) 



DCE in linear resonators. For linear cavities the 
E operator for transverse modes is 



E==*y]Cj(aj-aI)V'j(a;)e 



(7) 



with / = w^ — 4n^n^ ^ u}q{1 — 2Sn h{t)/no), by intro- 
ducing the variable y = E e^" and considering small 
modulations of the refractive index (Sn/riQ ^1). 

For periodic modulations / with period T, the formal 
solution of eq. (pi) is 



m 



S{t) 



y{t = 0) 
i/it = 0) 



(3) 



According to Floquet's theory [13], to characterise the 
evolution of the system at times t = mT (with m an 
integer) it suffices to consider the Lyapunov exponent fi 
(i.e. the rate of growth of the amplitude of the solution) 
which is computed as 



p = T-^ cosh-\Tr[2S{T)]}. 



(4) 



The solutions are bounded if \Tr[2S{T)]\ < 1, and un- 
bounded for |Tr[2S'(T)]| > 1, as a result of a parametric 
resonance (PR). Including optical losses, the photon yield 
at a PR is 



N{mT) 



2eono 

hujo 



m 



,2(m-7)™T 



(5) 



where E^ is the initial field amplitude and 7 is the coup- 
ling of the cavity with the environment. 

At T = 0, the amount of light emitted from the ground 
state of the field depends on the nature of the field con- 
sidered: in the classical model (-Eq) — and there is no 
photon emission whereas in the quantum model (-Eq) ^ 
due to the ZPF, yielding an exponential amplification of 
the field. But for T > 0, the difference between the clas- 
sical and quantum prediction is very small because the 
contribution from thermal fluctuations present in both 
models is relevant and masks the contribution from ZPF. 
At T > 0, the quantum character of the DCE is hid- 
den in the quantum correlations, calculated using a fully 
quantum dynamics. This is done replacing the E field 



where x is the coordinate along the optical axis of the res- 
onator, ^i{t) — •\/Sfcoc|j|/2rt2(i), a and Oj are the cre- 
ation and annihilation field operators (with [oj, a/-] — 5ij), 
Gj is the polarisation unit vector and tpjix) — sin{kjx) 
with kj — TTj/L (with j a positive integer). 

Replacing ^ into (fTl) yields the evolution eqs. for 
a, and Uj, which for each mode ipj correspond to the 



Hamiltonian 



H{t) - /(i)a]a, 



W (t) [af - o. 



(8) 



with f(t) = l/n{t) and g{t) = i/2(i In n(t)/dt. From this 

point all eqs. are written in natural units of the vacuum 

frequency, defined by h : 

For Gaussian states 



1 and kjC = 1. 



the Wigner function has 



the form W{X) 



(27r) det [a] ^ exp{- 



i(X 



X) a ^{X — X)}, where X = [x,p] and ct is a co- 
variance matrix defined by cr = XX"^ — X X^ . For the 
Hamiltonian ((sl), the master eq. ^ yields the evolution 
eqs. 



X = -(Mj+7/2)X, 
(7 = —M2 (7 — (JM2 — ^ [a — (Ja 



where 



M, = 2 



-.9 / 
-/ 5 



(9) 
(10) 



(11) 



and where 0^0 = (2?t. -f l)Id2/4 and Id,i is the n- 
dimensional identity matrix. The state cr 00 describes the 
cavity in equilibrium with the bath in the absence of the 
refractive index modulation. Notice that CTqo is de facto 
the ground state of the cavity at T > therefore, it is 
simultaneously the initial state of the cavity at t = (i.e. 
"■q = o'oo) and the asymptotic state to which the cavity 
decays if the modulation if turned off. 



DCE in ring resonators. For ring cavities the E 
operator for transverse modes is 



E ^^Y.^At) \aAt)4>j{^) - «]W0K^) 



L.J, 



(12) 



where (pji^) = exp{ikjx), with k±j = ±2Trj/L and j an 
integer (+j for hght propagation along the optical axis 
and —j in the opposite direction). In this case, there is 
a coupling between modes with symmetric wave number 
j resulting in the Hamiltonian [inj 

Hit) - f{t)[a]aj+al^a^,]+igit) [ajal^^a.a^,], (13) 
Making X' = TX, with X = [xj,pj,x-j,p-j] and 
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the evolution eqs. for X and a become 



X' = -{Mi 



-7/2)X', 



(14) 



(15) 



Ml a' - cr'Af4 - 7/2 [aa' + a' a - 2acr;„)(,16) 



Figure 2: Phase resonance 
The shaded regions are those where there is net amplification of 

the field (black: /,. = 1.2, grey: /^ = 1.5). Maximum 
amplification is obtained for Qi = ^/2 + kn but slight deviations 

from these values still produce amplification. The allowed 

deviation increases with the relative variation of the refractive 

index. 



where M4, a and a'^ are 4x4 matrices with a'^ ~ 
(2n + l)a74 and 



M, 



M2 








, a 



Id2 




(17) 



Notice that M4, a and a'^ are block diagonal. Then, 
eqs. (15 1 and (16 1 imply that, in terms of quadratures 
X' , two coupled modes of ring resonator can be decom- 
posed into two independent modes of a linear resonator. 
Moreover, one these modes is isolated from the bath and 
constitutes a DFS. Hence, when a PR occurs, it grows ex- 
ponentially with /i, regardless of the losses. This feature 
of the ring geometry does not exist for linear cavities, 
where all modes are coupled to the bath. 

Evolution of the state of the system. The solu- 



tions of eqs. ( 15 ) and ([l6|) for a ring resonator are 



X'{t) = UthX'o, 
a'it) = Uth<7'oUl + > 



(18) 
(19) 



where Cg = a'^ and Xq describe the initial state of the 
cavity, Uth = e~"'/^"*[/, U solves the eqs. in the absence 
of the bath and cr' is a particular solution of eq. ( 19 1. In 



this case, U and a' are 



and (T+ = 7C/+ Jq 



Mu^y 



dr 



ul 



. The matrix S is 
([2]). The solutions 



again obtained as the solution of eq. 

of eqs. ([9]) and ( 10 1 for a linear resonator are formally 

identical to those of eqs. ( 15 ) and ([161) '^'-"^ ^^^ degree of 



freedom coupled to the bath expressed in terms of X' . 

Examples of refractive index modulations. The 

impact of the DFS is clearer by considering the two par- 
ticular optical modulations: (i) a periodic and instantan- 
eous change of n between two values, rii and n^, which is 
directly related to TR; and (ii) a sinusoidal modulation, 
similar to more conventional setups of DCE [8 . 

For a sequence of instantaneous perturbation of the 
refractive index with period T given by 



fit) 



/i, rriT < t < mT + <i, 

/2, mT + ti<t<mT + ti+t2 = {m+l)T, 

(21) 



Using eqs. (15) and (16 1, it results that 



Tr[2S{T)] = 



(1 



4/,. 



cos 6'+ —„ cos ( 



4./r 



(22) 



U 



U+ 

u_ 



cr+ 




where U± = e 
/(i = 0). 



i?+ = 



±GS D-1 



(20) 



R±^SR± withG = !l9{T)dT, /o 



r 1] 




r 1 1 




r 1 1 


-/o 


,i?_ = 


0/0 


,s = 


-1 



with 0± - 01 ± 02, Oi = fiti, 02 = /2t2 and fr = /2//1. 
For 7 = 0, eq. (([22])) shows that the highest photon 
emission rate occurs for 0i = ^/2 + kir (with k integer) 
and fi = T~^ |ln/r|. For the same 9i but with 7^0, 
the asymptotic values of the number of photons emitted 
Af +1 = J2i=i 2i'H'^i + h.c.)/2 and of the entanglement 
of the state (measured using logarithmic negativity) at 



t = mT — > oo are respectively 
2n+l 



(JV + l) 

EN{mT) 



g2pmT_^g-2,,_mT 



mT-i-oo 4 

+Fl(m)+F;(m)] 

/iTTlT 

— > maxjO, - — 

mT— >oo In 2 



(23) 



logs 



{21% + 1) FUmy/A} (24) 



o27«l 



s27*ifl 



with F4(m) = [1 - e-2';±™T][i 
e2^*^)/±i]/[l - e^"*^] and 77± = (7 ± //). 

In the case of a small sinusoidal perturbation of the 
refractive index, given by n = no + Snsinnt, and for 
6n/no ^ 1, yields that /^ ~ [1 — 2((5n/no) sin nt]/nQ and 
eq. (pi) reduces to Mathieu's eq. 

j/+((5 + esin2i)?/ = (25) 

with t' = m/2, S = i/nln^ and e = -85n/nln^. Eq. 
( 25 1 exhibits PR for S = m^ (with m integer) or altern- 
atively fl = 2/mnQ, yielding ^, — Sn/2nQ. 

The asymptotic values of M and Ej\j at t = mT — )• cxi 
are 

(AA + 1) ^ ?^^ [e^'^™^ + e-2'(-™^ + 

mT— i-oo 4 



En 



+Fl{m) + F^{m)] 
^.mT 



(26) 



— > max{0, 

mT— >-oo In 2 



logs 



[2n + l)FUmy/A} (27) 



±(to) = [27(1 - e2''TT)(i 



'2r;±mT 



)]/Nt(1 



with F: 
e^^i^)]". 

In both cases there are two distinct contributions 
for M. The terms in (231 and (26) proportional to 
e^''™"^ result from the degree of freedom associated with 
the DFS, resulting in a photon emission rate of tdfs = 
2/i, independently of the bath or losses. The terms pro- 
portional to e^^''-™^ and F_, result from the remaining 
degree of freedom which is coupled to the bath and has 
two regimes: (i) a strong losses regime for rj- > 0, when 
there is no exponential amplification of the field and (ii) 
a weak losses regime for 77- < 0, when photons are pro- 
duced at the rate Tth = 2(fj, — 7). This is the emission 
rate for linear resonators, in which there is no DFS. 

The other relevant feature of the solution is the ex- 
istence of entanglement in both regimes. In fact, since 
i^-f (00) is constant, the linear growth term always dom- 
inates for large m. This means that a long enough mod- 
ulation of the medium will always result in entanglement 
after an occurrence time 



log2 [(2r7;+l)i^+(cx))V2]ln2 



(28) 



It is important to note that a finite temperature does 
not forbid the existence of entanglement but delays its 
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Figure 3: Occurrence time 

The occurrence time as a function of the temperature for different 

values of 7. The occurrence time increases with the temperature 

and with 7. 



occurrence, which enforces the need to work with low 
temperatures. These results compare well with reference 
|13) . where the author considered two harmonic oscil- 
lators prepared in an initial two mode squeezed state, 
not directly coupled but sharing a thermal bath. It was 
shown that, if the initial state was sufficiently entangled, 
it would preserve some amount of entanglement after any 
arbitrarily large time. In the case presented here the ini- 
tial state is disentangled and it is the DCE that produces 
correlations until they reach a threshold imposed by the 
thermal bath, after which entanglement is present. The 
common factor in both cases is the existence of a DFS 
which shields the entanglement from the environment. 

Conclusions. This letter introduces the idea of using 
entanglement as a marker of the quantum nature of the 
DCE at finite temperature and proposes a new experi- 
mental setting which has a DFS and is capable of sup- 
porting and amplifying entanglement. In more detail, it 
was shown that for a resonant modulation of the refract- 
ive index of the dielectric medium in a ring resonator it 
is always possible to amplify radiation, regardless of the 
loss regime. Furthermore, it was verified that the exist- 
ence of a decoherence free subspace allows quantum cor- 
relations to survive, making entanglement viable even at 
high temperatures. Entanglement is a distinctive charac- 
ter of a quantum effect which can be used to sort out the 
DCE from classical parametric amplification of thermal 
fluctuations. Hopefully these results will pave way to the 
observation of the DCE and entanglement at finite tem- 
peratures. 
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